To Appear in ApJ 

Preprint typeset using I^T^jX style cmulatcapj v. 10/09/06 

SPITZER 24/xm TIME-SERIES OBSERVATIONS OF THE ECLIPSING M-DWARF BINARY GU BOOTIS 

Kaspar von Braun 1 , Gerard T. van Belle 1-4 , David R. Ciardi 1 , Mercedes Lopez-Morales 2 , D.W. Hoard 3 , and 

Stefanie Wachter 3 

To Appear in ApJ 

ABSTRACT 

We present a set of Spitzer 24/um MIPS time series observations of the M-dwarf eclipsing binary 
star GU Bootis. Our data cover three secondary eclipses of the system: two consecutive events 
and an additional eclipse six weeks later. The study's main purpose is the long wavelength (and 
thus limb darkening-independent) characterization of GU Boo's light curve, allowing for independent 
verification of the results of previous optical studies. Our results confirm previously obtained system 
parameters. We further compare GU Boo's measured 24/im flux density to the value predicted by 
spectral fitting and find no evidence for circumstellar dust. In addition to GU Boo, we characterize 
(and show examples of) light curves of other objects in the field of view. Analysis of these light curves 
serves to characterize the photometric stability and repeatability of Spitzer 's MIPS 24/im array over 
short (days) and long (weeks) timescales at flux densities between approximately 300-2, 000/jJy. We 
find that the light curve root mean square about the median level falls into the 1-4% range for flux 
densities higher than lmJy. Finally, we comment on the fluctuations of the 24/im background on 
short and long timescales. 

Subject headings: techniques: photometric, infrared: stars, stars: fundamental parameters, stars: 
binaries: eclipsing, stars: individual: GU Boo, circumstellar matter, dust 



1. INTRODUCTION 

GU Bootis is a nearby, low-mass eclipsing binary 
system, consisting of two nea rly equal mass M-dwarfs 
([Lopez-Morales fc Ribasl 12005( 1 . It is one of currently 
very few (~5) known nearby (< 200 pc) double-lined, 
detached eclipsing binary (DEB) system s composed of 
two low- mass stars ((Lopez-M orales 2007). Eclipsing bi- 
naries can be used as tools to constrain fundamental stel- 
lar properties such as mass, radius, and effective tem- 
perature. Given the fact that over 70% of the stars in 
the Milky Way ar e low-mass objects with M < IMq 
(|Henrv et al.lll997l ). coupled with the considerable uncer- 
tainty over the mass-radius relation for low-mass stars, 
objects such as GU Boo are of particular interest in ex- 
ploring the low-mass end of the Hertzsprung-Russcll di- 
agram. 

While simultaneous analysis of DEB light curves and 
radial velocity (RV) curves provides insight into the 
component masses and physical sizes, an estimate of 
their intrinsic luminosities can only be made with the 
knowledge of the amount and properties of dust along 
the line of sight. As such, it is important to un- 
derstand whether low-mass DEB systems used to con- 
strain stellar models contain dust which, in turn, may 
lead to an underestimate of their surface temperatures 
and thus l uminosities. This pr o blem has bee n docu - 
mented in iDelfosse et all ((1999( 1: iMazeh et all (|2001( 1: 
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[Torres fe Ribasl (|2002l ): iRibad ((2003( 1. In particular, 
I Ribasl (|2003( l states that the most likely explanation for 
the temperature discrepancy between observations and 
models for the low-mass DEB CU Cancri is the presence 
of either circumstellar or circumbinary dust. The detec- 
tion of dust in any system such as GU Boo would there- 
fore additionally shed insight into formation and evolu- 
tion of the low-mass DEBs. 

The characterization of the effects of limb darkening 
and star spots introduces additional free parameters and 
thus statistical uncertainty in the calculation of the stel- 
lar radii and masses. Using the Spitzer Space Telescope, 
we obtained 24/im time series observations of three sepa- 
rate instances of GU Boo's secondary eclipse (see 32. 1[) to 
create a light curve far enough in the infrared to not be 
contaminated by the effects of limb darkening and star 
spots. We purposely timed the observations such that 
each secondary eclipse event is preceded by a sufficient 
length of time to establish GU Boo's out-of-eclipse flux 
density in order to detect any infrared excess possibly 
caused by thermal dust emission. 

A further goal of our study is to characterize the 
photometric stability of the Multiband Imaging Pho- 
tometer (MIPS) on Spitzer at 24/im over short and 
long time scales, similar to what w as done for bright 
objects in §5 of iRieke et all ((2004( 1. Time -series ob- 
serving is atypical (albeit increasingly common) for 
Spitzer, which is the reason why there are very few pub- 
lished photometric light curves based on Spitzer obser- 
vations. The recent spectacular observations of primary 
and secondary eclipses of tra nsiting planets are notable 
exceptions (see for i nstance [C harbonncau ct al. 120051: 
I Deming et al.1 120051: ICowan et all l2007t IDeming et alJ 
l2007t iGillon et all 120071 .iKnutson et al.ll2007D . Of these, 
the IDeming et al.l ((20051 ) study was performed at 24/im. 
We therefore observed two consecutive secondary eclipses 
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of GU Boo (~ 12 hours apart), and then a third event 
about six weeks later (see Tabled]). 

We describe our observations and data reduction meth- 
ods in ij2] and discuss our findings with respect to 
Spitzer 's photometric stability in $3l The analysis of GU 
Boo's light curve is described in $4] We probe for the ex- 
istence of an infrared excess in GU Boo's spectral energy 
distribution in In §6l we show light curves of other 
well sampled objects in the field along with a brief sum- 
mary of their respective properties, and we summarize 
and conclude in [J7] 

2. SPITZER OBSERVATIONS AND DATA REDUCTION 
2.1. Observations 

We used the MIPS 24/xm array aboard the Spitzer 
Space Telescope (|Werner et al.ll2004r i to observe GU Boo 
in February and April of 2006, as outlined in Table [TJ 
The MIPS 24^m array (MIPS-24), is a Si: As detector 
with 128 x 128 pixels, an image sca le of 2.55" pixel" 1 , 
and a field of view of 5.4' x 5.4' (|Rieke et all |2004[ ). 
Our exposures were obtained using the standard MIPS 
24/im small field photometry pattern, which consists of 
four cardinal dither positions located approximately in a 
square with 2 arcmin on a side. For two of these cardi- 
nal positions, there are four smaller subposition dithers 
(offset by ~ 10 arcsec), and for the other two cardinal 
positions, there are three such subpositions. This results 
in a dither pattern in which Spitzer places the star at 14 
different positions on the array. 

Our goal was to observe three independent secondary 
eclipses of GU Boo: two consecutive ones and another 
one s everal weeks after the first two (|von Braun et al.l 
2007). Of our total of nine of Spitzer 's Astronomical 
Observation Requests (AORs), three were used for each 
secondary eclipse event (see Table [lj . Each AOR con- 
tained eight exposures 5 with 36 individual Basic Cali- 
brated Data (BCD) frames each. The first BCD in each 
exposure is 9s long, the subsequent 35 are 10s long. The 
first two BCDs of every exposure were discarded due to 
a "first frames effect" . This procedure left 34 BCDs per 
exposure, 272 BCDs per AOR, 816 BCDs per secondary 
eclipse event, and 2448 BCDs for the entire project (all 
10s exposure time). 

For background information on Spitzer 
and MIPS, we refer the reader to the 

Spitzer Observer's Manual, obtainable at 

http : //ssc . spitzer . caltech . edu/documents/som/. 
For information specifically related to MIPS data reduc- 
tion, please consult the MIPS Data Handbook (MDH 
- http: //ssc .spi tzer . caltech . edu/mips/dh/) and 
IGordon et, aD ^K)M ). 

2.2. Data Processing and mopex / apex Reduction 

2.2.1. Mosaicing 

The MIPS-24 data are provided by the Spitzer Archive 
in the (flatfielded) BCD format. We applied further post- 
processing to these data in order to correct for small scale 
artifacts, in particular using IRAF's 6 CCDRED package 

5 The term "exposure" here can be thought of as a cycle of 
observations, but since the word "cycle" is reserved for another 
unit of Spitzer data collection, it carries the somewhat misleading 
name "exposure". 

6 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 



TABLE 1 

Spitzer. MIPS-24 observations of GU Bootis 



Date (zUUoJ 


A T T T") C f ' : 

M1FS Campaign 


Obs. Set AUHs 


Exposures a 


Feb 20 


MIPS006500 


1 16105472 


860 






16105216 








16104960 




Feb 21 


MIPS006500 


2 161041 04 


860 






16104448 








16104192 




Apr 01 


MIPS006700 


3 16103936 


860 






16103680 








16103424 





Note. — Two consecutive secondary eclipses were observed in observing 
sets 1 and 2, and a third secondary eclipse event six weeks later in observing 
set 3. 

a 10 seconds per exposure. 




Fig. 1. — A Spitzer MIPS 24/xm mosaic of GU Boo (marked 
with a circle at the center of the image). This mosaic was created 
using all 272 frames in one AOR and is about 8 arcmin on a side. 
North is up, east is to the left. The change in noise structure as 
a function of position is due to different effective exposure times 
(only the inner ~ 3' X 3' were covered by all 272 BCD frames). The 
white specks in the NW corner are flatfielding residuals, caused by 
a fleck of paint or dust grain on the pickoff mirror (from Spitz er 's 
launch), imaged at the four cardinal dither positions (see H2.lt . 

to remove the weak "jailbar" features in the images (as 
described in the MDH). 

The Spitzer software package m opex (|Makovoz fc Khar] 
[20051: iMakovoz k Maxleaull2005h was used to co-add the 
individual MIPS BCD frames into mosaics of 17 frames, 
using overlap correction and outlier rejection in the pro- 
cess. The choice of 17 frames was made to balance three 
aspects: 

1. We want to obtain sufficient signal-to-noise ratio 
(SNR) for measured flux densities in the combined 
images and subsequent data points in the light 
curves (SNR > 10 for GU Boo; see Table©. 

2. We need to maintain a sufficiently high effective 
observing cadence to temporally resolve elements 
of GU Boo's light curve for fitting purposes. 

3. We do not want to be forced to combine frames 
from different exposures into a single light curve 
data point (see §2.ip . 

Research in Astronomy, Inc, under cooperative agreement with the 
National Science Foundation. 
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The interpolated, remapped mosaics have a pixel scale 
of 2.45" pixel -1 . We show in Figure [l] the MIPS-24 field 
of view of GU Boo. 

2.2.2. Photometry 

For photometric reductions of the mosaiced images, we 
utilized the apex component of mope x to perform point- 
source extraction as described in iMakovoz fc Marleaul 
(2005) 7 . This step included background subtraction of 
the images, and the fitting of a resampled point response 
function (PRF). In order to match the PRF centroid as 
closely as possible to the centroid of the stellar profile, 
the first Airy ring is initially subtracted from the stellar 
profiles, and the source detection happens on the result- 
ing image. Photometry of the detected sources is then 
performed on the original images. 

For single frame photometry on mosaiced images, apex 
provides the option of using a template model PRF pro- 
duced by the analysis of ~ 20 bright, isolated stars in 
the Spitzer Archive, or using one's own data to create 
a PRF. It furthermore allows for a resampling factor in 
both x and y directions. The scatter in our light curves 
was minimized when using the model PRF provided by 
the Spitzer Science Center, oversampled by a factor of 4 
in both x and y directions. Using the PRF created from 
our own data or sampling any PRF to a higher resolution 
resulted in noticeably larger root mean square (rms) dis- 
persion in our light curves, most likely due to systematic 
errors introduced in the low SNR regime of our data (see 
Table O. 

We note that, currently, apex only provide s the 
optio n of using a synthetic PRF (Tiny Tim 8 ; iKristl 
1993) for photometry on ind ividual BCD f r ames (i.e., 
not mosaiced), as appli ed by iDeming et alj (|2005l ) and 
iRichardson et al.1 (|2006f) . Since the flux density of our 
target sta r (~ 600/zJy) is so m uch lower than HD 209458 
(~22mJv; foeming et aTll2005ft . our SNR regime did not 
allow for performing photometry on single BCD frames. 

2.2.3. Background Fluctuations 

The apex error analysis is described in the apex 
User's Guide, and pa rts of it can be found i n 
Makovoz et all (l2002ft j__ jMakovoz fc Lowrancd (|2005h : 



Makovoz fc Marleaul (|2005). We briefly summarize the 



general idea here. Errors in the photometry are domi- 
nated by the statistical background fluctuations in the 
images. These fluctuations are calculated per pixel by 
estimating the Gaussian noise inside a sliding window 
whose size is defined by the user (45 x 45 interpolated 
pixels in this case). Thus, apex produces "noise tiles" 
for the computation of the SNR of the point sources in 
the corresponding mosaiced image tiles (see column 7 in 
Table EJ. 

To provide an estimate of the background fluctuations 
from image to image, we show in Fig. [3 the surface 
brightness for every image in the three observing sets. 
These estimates were obtained by calculating the me- 
dian surface brightness level for the inner 90% of the 
image (in area). The error bars correspond to the stan- 
dard deviation about this median over the same area. 

7 Also see information on apex at 

http://ssc.spitzer.caltech.edu/postbcd/apex.html and the 



Note that the fluctuations of the background within ob- 
serving sets are very small, but they are different for the 
temporally offset observing set 3 (see Table [1]). Surface 
brightness values are given in the Spitzer native units of 
MJy per steradian. The surface brightness values in the 
cores of the brightest objects in the field are typically 1- 
1.5 MJy/sr above the background level (23-25 MJy/sr). 
A linear fit (weighted by the standard deviation values 
of the data points) to observing sets 1 & 2 returns a 
slope of -0.033 ± 0.023 MJy/sr/day. The same fit for 
all three observing sets produced a statistically consis- 
tent slope of —0.030 ± 0.001 MJy/sr/day, indicating a 
smoothly decreasing background level over the course of 
our observations (Table [T|) . 

The Spitzer tool Spot 9 predicts the surface bright- 
ness of the Zod i acal background in Spitzer images. 
lHarrington et al.1 (|2006h . for instance, used this model 
to correct for ostensible variations in instrument sensi- 
tivity over a period of around four days. We compared 
our background estimates to the predictions in Spot and 
found that the model underestimates our measurements 
by 3.5-4.0 MJy/sr. One potential reason for any off- 
set between observed and estimated backgrounds is that 
Spot calculates a monochromatic background, whereas 
the measured background is integrated over the wave- 
length passband and sensitivity function of the MIPS-24 
detector. Using Spot, we calculated background esti- 
mates in 6 hour increments from 2006 Feb 20 00:00:00 
to 2006 Feb 22 00:00:00, and from 2006 Apr 01 00:00:00 
to 2006 Apr 02 00:00:00. The slopes of the Feb back- 
grounds and the Feb-Apr backgrounds are -0.022 and 
-0.032 MJy/sr/day, respectively; Spot does not provide 
error estimates in its predictions. It thus appears that, 
within statistical uncertainties, the behavior of the Spot 
background model is consistent with our empirical results 
(at least for the time scales of our obse rvations) , lending 
furthe r justification to the approach bv lHarrington et alJ 
(|2006l) . 

If, however, the difference in slope between the Spot 
estimates for just Feb and Feb-Apr is indeed real but 
simply not detectable at our temporal resolution (in- 
dicating that the background change with time is not 
a simple linear function), any resulting discrepancy be- 
tween calculated and observed background values may 
be attributable to the fact that the model is calculated 
for Earth, whereas Spitzer is in an Earth-trailing orbit 
and t hus looking through di fferent amounts of Zodiacal 
dust ([Harrington et al 



i gn di rt 
, 2006). 



3. PRECISION AND REPEATABILITY OF THE SPITZER 
PHOTOMETRY 

3.1. PRF Fitting Versus Aperture Photometry 

To create GU Boo's 24/im light curve, we performed 
both PRF photometry as described in £12.21 and addi- 
tionally utilized apex's option of simultaneously obtain- 
ing aperture photometry. Figure [3] shows the agree- 
ment between PRF and aperture photometry based on 
an aperture radius of 6" (with 20-32" background annu- 
lus) which minimized the rms in the flat part of GU Boo's 
phased 10 light curve. Using information from the MDH, 



User's Guide at http: //ssc . spitzer . caltech. edu/postbed/ doc/apex .pdf 



9 http: //ssc . spitzer . caltech. edu/propkit/spot/index .html 

l 1 " We use the period calculated by Lopez-Morales & Ribas 



3 http: //ssc . spitzer . caltech. edu/archanaly/contributed/stinyti|i2.005) (0.488728 days; see Table [2]l for our phasing throughout 
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Fig. 2. — Temporal fluctuations of the image background mea- 
sured in surface brightness units of MJy/sr for the three observing 
sets JTJl. The lower, middle, and upper panel correspond to observ- 
ing set 1, 2, and 3, respectively. Note that the background level 
is similar between observing sets 1 and 2 (24.48 ± 0.009 MJy/sr 
and 24.46 ± 0.009 MJy/sr, respectively), but lower for observing 
set 3 (23.28 ± 0.008 MJy/sr). The typical fluctuations within a 
given image (represented by the error bars) are around 0.05-0.06 
MJy/sr. The large error bar for the one data point in observing 
set 3 is caused by a cosmic ray. 

we applied a multiplicative aperture correction of 1.699 
to the photometry. The median flux density obtained 
by PRF photometry for the flat part of the phased light 
curve is 614 ± 49/iJy compared to 608 ± 59/iJy for the 
aperture photometry. Thus, the absolute median flux 
density values agree very well for the two different pho- 
tometry approaches, but the PRF photometry exhibits 
smaller scatter around the median magnitude. We note 
that the current version of apex does not calculate pho- 
tometry errors in the aperture correction, and the prin- 
cipal reason why we performed aperture correction is to 
verify the absolute flux density level of our sources as 
calculated by PRF fitting. 

3.2. Absolute Versus Relative Photometry 

In order to remove statistically correlated noise from 
GU Boo's light curve, we performed re lative photometry 
as des cribed in equations 2 and 3 of lEverett fc Howelll 
HQOH). We picked comparison objects based on the num- 
ber of observational epochs: in order to obtain a rela- 
tive offset per photometric data point in GU Boo's light 
curve (all data points are treated independently of each 
other), it is advantageous to use stars with (at least) 
as many data points as GU Boo itself. Four objects 
out of Table [3] fulfill this criterion: numbers 18, 31, 58, 
and 66 (see Figures [TOHTSl for their light curves) . The 
cross-referencing in Table |3] shows that objects 31 and 
66 are stars, and objects 18 and 58 are galaxies (as are 
all other objects in the field that we were able to cross- 
reference VJNotejJrowevei^^ 31's closest match 
in SPSS lAdelman-McCarthv & et alll200l and 2MASS 
(|Cutri et alj|2003al ISkrutskie et al.ll2006D catalogs is 9" 
away whereas for object 66, the distance to the closest 
SDSS match was only 0.15". We originally presumed 
that, despite the fact that they are galaxies, objects 18 
and 58 would be unresolved at our large pixel size ( §2.ip 

this paper. 
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Fig. 3. — Comparison between PRF photometry (blue crosses) 
and aperture photometry (red squares) for GU Boo's phased light 
curve. Photometric error bars are omitted for the sake of clar- 
ity. The scatter in the flat part of the light curve is smaller for 
PRF photometry (49/iJy) than aperture photometry (59^tJy), but 
the median flux density level of the flat part of the light curve is 
identical within the errors for both photometry approaches. 















□ 




□ 


□ 






















□ 






, i , , , , i 


, , 1 , 


, t i . . i T . i 





0.1 
phase 



0.1 
phase 



Fig. 4. — Comparison between PRF photometry (blue crosses) 
and aperture photometry (red squares) for objects 31 (left) and 66 
(right). Photometric error bars are omitted for the sake of clarity. 
The light curves are phased to the period of GU Boo for the sake 
of com parison. The ordinate scale is the same as Figures [3] |5] and 
I10H12I The flux density values (PRF vs aperture photometry) are 
identical within the errors for both objects. 

and tested that hypothesis by comparing the flux den- 
sity obtained by PRF photometry to that obtained by 
aperture photometry (6" aperture; ^3.ip . Figures 2] and 
[5] shows that our presumption does not hold true for ob- 
ject 58, and we discarded it from our relative photometry 
procedure. 

We find that, by performing differential photometry as 
outlined above, the scatter in the flat part of GU Boo's 
phased light curve reduces by 9.4% over the PRF pho- 
tometry (see Fig. \S§ to 45/zJy. Light curve fitting as 
described in ||4] was performed on the differential pho- 
tometry. 

3.3. Intra-Set Versus Inter-Set Photometric Stability 
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Fig. 5. — Comparison between PRF photometry (blue crosses) 
and aperture photometry (red squares) for objects 18 (left) and 58 
(right). Photometric error bars are omitted for the sake of clarity. 
The light curves are phased to the period of GU Boo for the sake 
of comparison. The ordinate scale is the same as Figures [3] [4] 
and !10lfT2l The flux density values (PRF vs aperture photometry) 
agree within the errors for object 18, but are discrepant for object 
58. 

Figure[6]shows fractional rms values versus median flux 
densities for all objects in Tabled For every object, we 
plot fractional rms for each individual observing set as 
well as for the three sets combined. Observing sets 1 
and 2 were obtained during the MIPS006500 campaign, 
observing set 3 during MIPS006700 (Table Consis- 
tent with the results in Rie ke et"aT] (|2004f ) . we find that 
inter-set repeatability of Spitzer 's MIPS-24 is compara- 
ble to the intra-set repeatability, both in terms of median 
flux density and the rms scatter of the light curves, de- 
spite varying background levels (see ij2.2.3p . For objects 
with a flux density in excess of lmJy, the rms scatter 
approaches 1-2 %, similar to the scatter fo und for the 
bright est sources observed with MIPS-24 in Riek e etldl 
(2004). Because of the intrinsic variability produced by 
the stellar eclipse, GU Boo has the largest fractional rms 
(~ 0.2). However, when we subtract the fit (see @]) from 
GU Boo's light curve, the fractional rms falls to 0.081, 
consistent with stars of similar median brightness. We 
show GU Boo's light curve for the three individual ob- 
serving sets in Fig. [7] and the phased light curve along 
with the fit in Fig. El 

In order to compare our rms values to background- 
limited noise values, we used Spitzer's SENS-PET 11 to 
predict the MIPS-24 sensitivity (la above background 
for 170 seconds integration time; see iJ2.2.1[) for low and 
medium background levels (solid and dashed line in Fig. 
O respectively). The mid-infrared background at the 
time of observations of GU Boo is 23-24.5 MJy/sr, which 
is between the typical low and medium background levels 
used by SENS-PET (see also N2.2.3p . We find that the 
SENS-PET predictions are consistent with our empiri- 
cally determined error estimates. Except for the vari- 
able GU Boo, typical values for the rms scatter of the 
light curves (Table [3]) are approximately equal to aver- 
age photometric measurement uncertainties of individual 
data points (see Figs. [THHU and [T2")). 

11 http: //ssc . spitzer . caltech. edu/tools/senspet 
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Fig. 6. — Median flux density versus fractional rms for the 24 
objects that have photometry for more than 72 out of 144 obser- 
vational epochs. Triangles, squares, and crosses represent the data 
from observings sets 1, 2, and 3, respectively (see Table [TJ, to 
illustrate the repeatability of S'pitzer/MIPS-24 within individual 
observing sets. Circles mark the data points from the combina- 
tion of all 3 observing sets (to show the inter-set stability). The 
data point with the highest fractional rms is GU Boo, due to its 
intrinsic variability. When subtracting our fit from its light curve 
(see Fig. [§}, GU Boo's fractional rms falls to 0.081. The solid and 
dashed lines indicate MIPS-24's sensitivity for our exposure times 
as a function of flux density for low and medium background levels, 
respectively. 
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Fig. 7. — GU Boo's 24/ mi lig ht curve, based on absolute PRF 
photometry (see EI3.1I and i|3.2H . The three panels represent the 
three MIPS-24 observing sets during which the individual sec- 
ondary eclipse events were observed. GU Boo is object number 
51 in our numbering system (see Table [3). 

4. ANALYSIS OF GU BOO'S PHOTOMETRIC LIGHT 
CURVE 

We modeled the secondary eclips e observations of 
GU B oo using the JKTEBOP code (jSouthworth et all 
I2004afbh. JKTEBOP is based on the original EBOP 
code (|Popper fc Etzell I"l981l : lEtzell Il981h . but with the 
additio n of the Leyenber g-Marquardt optimization algo- 
rithm ([Press et al.l I1992T ) to find the best fitting model, 
and also the implementation of a Monte Carlo simulation 
algorithm to determine robust uncer tainties in the fitted 
parameters ()Southworth et al]|2005h . 

The orbital period and initial epoch of the primary 
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eclipse wer e set to the values given in the ephemeris 
derived bv lLopez-Morales fc Ribasl (|2005h . We further 
fixed the mass ratio and the radius ratio of the stars, as 
well as the eccentricity of the system (e=0) to the values 
obtained in that work. We assumed no limb darkening ef- 
fects in the light cur yes, as expected for observations this 
far in t o the infrared (IClaret et al . 1995; Ri chardson et al.l 
[20061 iCiardi et all |2007t |Snellen| |2007t and references 
therein), and no significant gravitational darkening or re- 
flection effects, based on the spherical shape of the stars 
and the similarity in effective temperatures. All these are 
reasonable assumptions, based on the results of the study 
of GU Boo at visible wavelengths, and they are, in fact, 
hard to test in detail, given the photometric precision of 
the Spitzer light curve at this flux density level. 

In the absence of primary eclipse observations, to cal- 
culate the luminosity ratio of the system, we place a fur- 
ther constraint to the fit by fixing the value of the surface 
brightness ratio of the stars to J=J2/J1=0.9795. This 
value, combined with the adopted radius ratio and the 
no limb darkening assumption, gives a luminosity ratio of 
of L2/L1 = 0.9697, which is consistent with the expected 
value for GU Boo at 24/xm. 

The parameters initially left free in the models were: 
(1) the fractional sum of the radii, i. e., (i?i + i?2)/a, 
where R\ and R2 are the component radii, and a is the 
orbital separation, computed from the stellar masses and 
the orbital period of the system, (2) the inclination of the 
orbit i, (3) the amount of third light L 3 , and (4) a phase 
offset parameter </> (to account for small errors in the 
ephemeris). 

Our best model solution is illustrated in Figure with 
a reduced \ 2 of 1.7, and a mean fractional error per data 
point of 9.5% (cf. Fig. [6]). Formal errors in the fitted pa- 
rameters were derived using the Monte Carlo algorithm 
implementation in JKTEBOP for a total of 1000 itera- 
tions. We obtain a radius for the secondary component 
of i? 2 = 0.66 ± 0.02 Rq. Our value of the orbital incli- 
nation is i = 89.3 ± 0.8 degrees. Both values are slightly 
larger than the ones obtained bv lLopez-Morales fc Ribaa 
(I2065I) at optical wavelengths, R 2 = 0.62 ± 0.02 R Q and 
i = 87.6 ± 0.2 degrees. The two secondary radius esti- 
mates agree to within random statistical errors (1.4cr). 
In the case of the inclination, our value is not as well 
constrained as in the optical, since we lack a full light 
curve that includes a primary eclipse. We show our esti- 
mates for GU Boo's system parameters in Table [2] For 
the tested third light contribution, we obtain a value of 
L 3 = -0.04 ± 0.07, consistent with L 3 = 0. 

Finally, we find a phase shift of = -0.014 ± 0.001. 
Thi s phase shift is 1.5 time s large r than expected from 
the lLopez-Morales fc Ribasl (|2005l) (\A<f>\ = 0.009), but 
can still be attributed to uncertainties in the origina l 
period estimation. The lLopez-Morales fc Ribasl (2005) 
observations were conducted in 2003, near JD=2452733. 
The number of elapsed periods inbetween those obser- 
vations and o ur Spitzer AORs is about 215 0. The 1- 
a error in the lLopez-Morales fc Ribas (20~05T) period es- 
timate is 2 x 10~ 6 days, which accumulates to 0.0043 
days, about 0.009 in phase, in 2150 periods. Thus, 
the disc repancy we find corresponds to about 1.5a 
from the lLopez-Morales fc~R ibas (2005) ephemeris pre- 
dictions. We estimate that this offset is based on normal 
statistiscal errors. An alternative explanation would be 



TABLE 2 

GU Boo System Parameters 



Parameter 



Value 



Orbital Period (days) a 

Orbital Eccentricity a 

Mass Ratio (M2/M1) a 

Combined out-of-eclipse 24 jim flux (^tJy) 

Radius of Secondary Component (Rq) 

Orbital Inclination i (degrees) 



0.488728 ± 0.000002 
(fixed) 
0.9832 ± 0.0069 
614 ± 49 
0.66 ± 0.02 (0.62 a ) 
89.3 ± 0.8 (87.6 a ) 



1 ILopez-MoTal cs & Ribai J2005D 



fc 400 





Fig. 8. — Our best fit overlaid on top of the phased 24/^m 
data of GU Boo derived from our relative photometry of the sys- 
tem. The fit has a reduced \ 2 = l-7± 0.1. The bottom panel 
shows the residuals around the fit. The calculated system param- 
eters, which agree well with t he results from the optical study in 
Lopez-Morales & Ribas ( 2005), arc shown in Table [2] 

that a third body orbiting the system could cause this 
shift, but since (1) we show in fj5] that GU Boo's flux is 
consistent with its modeled spectral energy distribution, 
and (2) we calculate the third light component to be L3 
= 0, any such claim would be unsubstantiable with our 
data. 

Equation [T] shows the updated ephemeris equation of 
GU Boo by combining the se ven minima in table 5 of 
lLopez-Morales fc Ribasl ([2005) with the thr ee new min- 
ima presented in this work. 



T(Min I) = HJD2452723.981327(1) + 0.4887247(8) • E. 

Uncertainty digits are given in parentheses. E rep- 
resents the number of elapsed periods since the initial 
epoch, T(Min I) the time of primary eclipse minimum. 

5. COMPARISON BETWEEN EXPECTED AND MEASURED 
24/xm FLUX DENSITY OF GU BOO 

In addition to the relative photometry of GU Boo (Fig. 
[5J, we also performed absolute photometry as reported 
in 33.11 The 24/mi flux density of 614 ± 49/zJy was de- 
termined from the median flux level outside of eclipse. 
To test the accuracy of the absolute flux density level, 
we show in this section a spectral energy distribution 
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(SED) model between 0.11 and 35/im, scaled to the op- 
tical and near-infrared (NIR) magnitudes of GU Boo 
(|L6pez-Morales fc Riba"sll2005[) . 

The GU Boo system components are two M stars of 
nearly identical mass, temperature and radius. For our 
SED model, we assumed both stellar components to be 
M1V stars with effective tem peratures of 3800 K (e.g., 
lLopez-Morales fc Ribasll2005l ). The model SED was con- 
structed from the M1V 0.11 - 2.5/xm optical-NIR tem- 
plates of iPicklei (|1998l ) a nd the Spitzer 5 — 35/^m In- 
frared Spectrograph (IRS; iHouck et al.ll2004f) spectra o f 
GL 229A, an M1V (3800 IP star ijCushing et al.l[200fih . 
To build the SED model (Fig. [9|), the M1V optical-NIR 
template was scaled to GU Boo's optical-NIR flu x dens i- 
ties based on table 1 in lLopez-Morales fc Ribasl (|2005f) . 
To connect the Spitzer IRS spectrum to the optical-NIR 
template, we fit a power law of the form F v oc v n (see 
dashed line in Fig. [9]) to the IRS spectrum. We found the 
best-fit exponent to the power law to be n = 1.9. The 
IRS spectrum and the power law (extrapolated to 2.4^tm) 
were then scaled to the red edge of the optical-NIR tem- 
plate. The slope of the power law was maintained to 
ensure a continuous transition between the optical-NIR 
template and the IRS spectrum (see Fig. [9]). Note that 
only the optical and NIR flux densities were used to scale 
the SED model; i.e., the scaling does not utilize the 24/im 
data point. 

The SED model predicts a mid-infrared flux density for 
GU Boo of i^(24^m) w 650/iJy. The measured 24/xm 
flux density of GU Boo (614 ± 49/iJy) is within la of 
the predicted flux density, agreeing remarkably well with 
the simple SED model presented here. We conclude that 
the stellar components are solely responsible for the mid- 
infrared emission of GU Boo. 

The few M- and K-dwarf DEB systems studied 
to date (GU Boo included) reveal that many of 
the binary components have larger radii (by 10- 
20%) and cooler effective temperatures (by 100 K to 
hundreds of K) t han predicted by stel l ar evo l ution- 
ary models (e.g., SErres & Ribasl 120021: [Ribasl 120031 : 
lLopez-Morales fc Ribasl l2005t iLopez-MoralesI I2007D . 
Magnetic activity and metallicity can a ccount for the ra- 
dius discrepancy ([Lopez-Mo rales 200^) and, in principle, 
also for the temperature discrepancy. An alternative ex- 
planation for the temperature discrepancy, however, is 
the presence of dusty material around the systems. The 
excellent agreement of our observed mid-infrared flux 
density with the model SED suggests that there is lit- 
tle, if any, (warm) circumstellar dust in GU Boo, likely 
ruling out circumstellar dust as a viable explanation for 
discrepancies with the stellar evolutionary models. 

6. LIGHT CURVES OF SELECTED OBJECTS IN THE 
FIELD OF GU BOO 

In this Section, we present a brief summary of selected 
other light curves in the field of GU Boo, along with basic 
determination of spectral types of the objects identified 
as stars (see §3.2|) . We limit our selection to the three 
objects that were used to perform the relative photome- 
try (see Ej3.2[) . Figures [TOHT!?! display these light curves. 
They are all on the same scale with different zeropoints. 
Parameters for all objects with at least 72 out of the 144 
epochs are listed in Table [31 We do not show light curves 
for the rest of the field objects since they can essentially 




A [ M m] 

Fig. 9.— The SED o f GU Boo, ba sed on M1V 0.11 - 2.5/im 
optical-NIR templates of Pickles ( 1998) and Spitzer 5 — 35^tm IRS 
spectra of GL 229A (spectral type M1V; T eff = 3800 K). The 
dashed line respresents the interpolation between template and 
spectra. For details, see [[5] 
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0.65 0.7 0.75 1.15 1.2 40.7 40.75 40.8 

HJD-2453706 

Fig. 10.— The MIPS-24 light curves of object 18 (a galaxy) in 
the field of GU Boo. The three panels represent the three MIPS-24 
observing sets. For parameters, see Table [3] 

be described as flat lines with some scatter around the 
median magnitude, which is characterized by the values 
in Table H 

Spectral typing for the two stars (objects 31 and 66) 
was attempted by means of SED fitting of photometry 
availab le in the literature: both objects ha ve Sloan DSS 
ugriz ([Adelman-McCarthv &; et all l2007| ) data points, 
and st ar 31 additionally has 2MASS JHK, (ICutri et all 
l2003allbT ) and Johnson RI (jMonet et all l2003f) magni- 
tudes available for it. SED fits were performed using 
the se dFit program discussed in §3.1 of Ivan Belle et all 
(|2007t ). The best SED match for star 31 is an M3III gi- 
ant (xleduced ~ 1-6), whereas star 66's SED was found 
to be consistent with an A2V dwarf {Xreduced ~ 0-9). 
Note that one assumption we make here is that the cross 
referencing for object 31 is correct, despite the large dis- 
tance from its closest matches in the SDSS and 2MASS 
catalogs (Table [3j. 

7. SUMMARY AND CONCLUSIONS 
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31 



HJD-24537B6 

Fig. 11— The MIPS-24 light curves of object 31 (possibly a 
star) in the field of GU Boo. The three panels represent the three 
MIPS-24 obser ving sets. For parameters, see Table [3] If the cross 
referencing in i|3.2l is correct, object 31 is an M3III giant (see text 
in JBJ. 






66 



1.15 1.2 
HJD-2453786 



Fig. 12.— The MIPS-24 light curves of object 66 (a star) in 
the field of GU Boo. The three panels represent the three MIPS- 
24 observing sets. For parameters, see Table \3\ Our SED fitting 
indicates this to be an A2V dwarf (see text in ^6j. 

We used MIPS-24 onboard the Spitzer Space Telescope 
to obtain time-series photometry of the M-dwarf DEB 
GU Boo. Our observations cover three secondary eclipse 
events, two consecutive ones and an additional event six 
weeks later. Analysis of the photometry shows that the 
flux density values for aperture photometry and PRF 
photometry agree, and that the PRF photometry pro- 
duces smaller scatter in the light curve. This scatter can 
be further reduced by performing relative photometry 
based on three comparison objects in the field. We find 
that the repeatability of MIPS-24 photometry is consis- 
tent over all temporal scales we sampled: within an ob- 
serving set and on time scales of 24 hours and six weeks. 



Our mid-IR analysis of GU Boo's light curve is less af- 
fected by stellar surface features than its optical counter- 
part. The results we produce show very good agreement 
with the previously obtained system parameters based 
on optical and near-IR work. A comparison between GU 
Boo's flux density and its model SED based on stellar 
templates and IRS spectra shows no IR excess, leading 
us to the conclusion that no warm circumstellar dust is 
present in the system. 

Finally, light curves of other objects in the field in- 
dicate that the photometric stability of Spitzer 's MIPS- 
24 is comparable over short (hours to days) and long 
(weeks) time scales, despite fluctuations in the image 
mid-IR background on time scales of weeks. 
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Basic Parameters of Objects in the Field of GU Boo 
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flllX'j ( //.Tv^ 


epochs 


2 


15 


:21 


:43 


.27 


33 


:52:53.00 


_ 


J152143.31+335252.6 


Galaxy 


11.2 


729 ± 69.19 


743 


± 65.53 


729 ± 82.25 


704 ± 42.61 


114 


3 


15 


:21 


:41 


.20 


33 


:53:10.99 


- 


J152141. 26+335310.1 


Galaxy 


9.0 


549 ± 61.34 


537 


± 64.70 


542 ± 57.65 


560 ± 55.56 


72 


7 


15 


:21 


:57 


.76 


33 


:52:27.17 


- 


J152157.77+335226.9 


Galaxy 


6.6 


414 ± 48.48 


400 


± 46.08 


424 ± 43.55 


418 ± 52.85 


129 


12 


15 


:21 


:51 


.50 


33 


:53:55.18 


- 


J152151.61+335355.1 


Galaxy 


13.9 


820 ± 45.50 


819 


± 38.65 


811 ± 51.52 


827 ± 38.69 


127 


16 


15 


:21 


:45 


.61 


33 


: 54:46. 99 


15214563+3354466 


J152145.64+335446.4 


Galaxy 


31.4 


1710 ± 66.06 


1720 


± 46.81 


1710 ± 63.68 


1630 ± 31.41 


120 


17 


15 


:22 


:00 


.15 


33 


: 54:00.34 


- 


- 


- 


7.0 


383 ± 43.07 


374 


± 48.01 


384 ± 37.40 


384 ± 40.65 


134 


18 a 


15 


:21 


:49 


.41 


33 


: 54:53. 36 


- 


J152149.41+335452.4 


Galaxy 


13.8 


727 ± 43.08 


727 


± 47.02 


731 ± 45.03 


725 ± 33.90 


144 


19 


15 


:21 


:40 


.34 


33 


:55:27.20 


15214031+3355259 


J152140.34+335526.0 


Galaxy 


39.9 


2120 ± 34.89 


2120 


± 33.17 


2140 ± 36.68 


2100 ± 26.63 


72 


29 


15 


:21 


:40 


.75 


33 


:55:54.64 


- 


- 


- 


7.9 


410 ± 39.22 


391 


± 31.21 


405 ±43.46 


430 ± 30.42 


72 


31 a 


15 


:21 


:52 


.83 


33 


:55:13.38 


15215249+3355046 c 


J152152.53+335504.9 C 


Star? 


6.6 


338 ± 33.17 


338 


± 30.70 


331 ± 26.36 


342 ± 38.59 


144 


37 


15 


:21 


:46 


.94 


33 


:55:49.89 


- 


- 


- 


12.1 


599 ± 48.87 


610 


± 40.35 


601 ± 53.10 


572 ± 33.25 


116 


11 


15 


:21 


:54 


.15 


33 


:55:44.33 


- 


J152154.21+335544.1 


Galaxy 


6.8 


336 ± 40.41 


313 


± 34.81 


316 ± 39.37 


357 ± 31.34 


136 


46 


15 


:21 


:52 


.71 


33 


:55:53.94 


- 


- 


- 


6.2 


300 ± 32.19 


298 


± 30.99 


299 ± 30.75 


300 ± 33.70 


129 


51 b 


15 


:21 


:54 


.80 


33 


:56:09.35 


15215482+3356088 


J152154.83+335608.9 


Star 


11.0 


559 ± 111.90 


520 ± 109.45 


554 ± 105.12 


579 ± 113.62 


143 


55 


15 


:21 


:44 


.87 


33 


:56:48.86 




J152144.89+335648.4 


Galaxy 


31.8 


1600 ± 39.13 


1570 


± 32.05 


1580 ± 30.08 


1630 ± 33.78 


72 


58 


15 


:21 


:48 


.90 


33 


:56:48.06 


15214889+3356478 


J152148.87+335647.3 


Galaxy 


16.7 


832 ± 38.17 


826 


± 35.79 


821 ± 35.08 


851 ± 32.40 


144 


65 


15 


:21 


:59 


.48 


33 


:56:20.17 








5.8 


276 ± 31.80 


280 


± 29.77 


277 ± 29.81 


259 ± 32.21 


115 


66 a 


15 


:21 


:56 


.32 


33 


:56:37.66 




J152156.31+335637.6 d 


Star 


9.0 


442 ± 34.76 


437 


± 37.28 


446 ±32.79 


441 ± 31.39 


144 


70 


15 


:21 


:58 


.89 


33 


:56:33.23 








5.6 


268 ± 24.85 


271 


± 23.22 


279 ± 24.02 


260 ± 20.10 


73 


74 


15 


:21 


:49 


.62 


33 


:57:20.07 




J152149.59+335718.8 


Galaxy 


5.2 


272 ± 24.46 


265 


± 27.86 


272 ± 18.05 


273 ± 25.82 


77 


84 


15 


:22 


:03 


.07 


33 


:57:24.75 


15220305+3357239 


J152203.05+335724.2 


Galaxy 


17.2 


908 ± 45.81 


896 


± 34.68 


903 ± 36.81 


969 ± 29.64 


120 


85 


15 


:21 


:57 


.43 


33 


:57:45.84 








5.9 


314 ± 36.37 


318 


± 29.56 


303 ± 32.52 


315 ± 40.68 


105 


90 


15 


:21 


:58 


.86 


33 


:59:14.34 




J152158.83+335915.2 


Galaxy 


11.8 


777 ± 58.61 


766 


± 54.21 


781 ± 64.90 


789 ± 48.39 


120 


91 


15 


:21 


:51 


.67 


3.3 


:59:58.95 








7.3 


516 ± 58.57 


496 


± 55.89 


517 ± 53.49 


563 ± 52.59 


9.3 



z 
td 

H 
% 



No te. — Para meters for all objects in the field of view of GU Boo with at least 72 epochs (half the total number). Listed arc source ID, position, 2MASS l lCutri ct al. 2003a: Skr utskie et al.l 
2006) and SDSS ( Adclman-McCarthv & ct al 2007) cross-referenced IDs, and SDS S type (if available). Flux densities are median values, rms denotes the scatter around the median. SNR 
(column 7) represents the average SNR per data point in the light curve (see §2T2) . Flux density in column 7 indicates the overall median value, whereas the flux densities fluxi, flux2, and 
flux3, (columns 9 — 11) show the values for the 1st, 2nd, and 3rd observing sets (Tablc[^, respectively. The last column is the total number of data points for the object. 

a Used as comparison object for differential photometry ( ^3.2t .^ GU Boo. c This is the closest match in both SDSS and 2MASS catalogs, with a distance of ~ 9", rendering the cross referencing 
somewhat uncertain. 1 ^ The distance between object 66 and the best SDSS match is around 0.15". 



